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ABSTRACT 

We present the results of an analysis of ultraviolet observations of T Tauri Stars (TTS). By 
analysing emission measures taken from the literature we derive rates of ionizing photons 
from the chromospheres of 5 classical TTS in the range ^ lO^'^-lO^'* photons s^^, although 
these values are subject to large uncertainties. We propose that the He Il/C IV line ratio can be 
used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted 
by TTS. By studying this line ratio in a much larger sample of objects we find evidence 
for an ionizing flux which does not decrease, and may even increase, as TTS evolve. This 
implies that a significant fraction of the ionizing flux from TTS is not powered by the accretion 
of disc material onto the central object, and we discuss the significance of this result and 
its implications for models of disc evolution. The presence of a significant ionizing flux in 
the later stages of circumstellar disc evolution provides an important new constraint on disc 
photoevaporation models. 
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1 INTRODUCTION 

It is now well established that at an age of ~ 10^ years, the major- 
ity of stars are surrounded by discs that are optically thick at optical 
and infrared wavele ngths ( Strom et al. 1989; Kenyon & Hartmann 
[199?; Haisch. Lada & Lada 2001). These discs are deduced from 
their sub-millimetre emission to be relatively massive (typically up 
to a few per cent of a solar mass; |Beckwith et al. 1990) and are 
thus widely held to be potential sites of planet formation. At an age 
of > 10^ years, however, most stars are no longer associated with 
such massive circumstellar discs, although low mass 'debris discs' 
(believed to consist of dust replenished by cometaiy impacts) have 
been detected through submillimetre and scattered light imaging 
in some early type stars of this age (Mannin gs & Sargenll[T997[ 
^yatt. Dent & Greaves 2003). An outstanding unsolved question 
is: by what mechanism do stars lose their primordial discs as they 
approach the main sequence? 

Some clues about the nature of the disc dispersal mecha- 
nism are provided by the distribution of infrared colours among T 
Tauri stars (Kenyon & Hartmann 1995). These colours fall rather 
neatly into two categories: those that are explicable in terms of 
a stellar photosphere plus optically thick circumstellar disc and 
those that are compatible with purely photospheric emission. (Note 
that these two categories tend to coincide with the spectroscopic 
designations of Classical and Weak Line T Tauri Star, hence- 
forth CTTS and WTTS.) As a number of authors have noted 
ISkrutsk ie etai]|l99 0'; 'Hartigan et al."l990l; iKenvon & HartmaniJ 
ri995; Armitage. Clarke & Tout 1999) , there is a striking lack of 
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transition objects with colours intermediate between these two cat- 
egories, suggesting that disc dispersal must occur on a time-scale 
that is short compared with the time the system spends as either a 
CTTS or a WTTS. 

Such two time-scale behaviour is however not readily pro- 
duced by some of the most obvious mechanisms for disc dis- 
persal, such as viscous evolution or magnetospheric clearing 
( Armitage et al. 1999). These instea d give rise to rou ghly power 
law declines in disc properties ( Hart mann et alj|T998h . so that the 
time required to decline from any state is always of the same 
order as the time spent in the foregoing state. Detailed models 
of disc photoevaporation were constructed by Hollenba ch et alJ 
(19,94. 2000) for the cases of both weak and strong stellar winds, 
but these models also fail to produce the required two-time- 
scale behaviour. Subsequent studies have improved on the de- 
tails of the photoevaporative winds, studying effects such as dust 
scattering ( Richling & Yorke 1997) and more detailed hydrody- 
namics (.Font et al. ,2004) . However the qualitative behaviour re- 
mains unchanged from that predicted b y the initial model s of 
iHoUenbach et alj <1994 . Very recent work <Adams et alj2004 has 
shown that far-ultraviolet (FUV) photons, which can dissociate H2 
without ionizing neutral hy drogen, can also driv e a significant disc 
wind. However the study of lAdams et alj j2003) focuses on an ex- 
ternal radiation field, and it is not yet clear if a central FUV source 
can produce a similar effect. The only model that has so far suc- 
ceeded in achieving the required two time-scale behaviour com- 
bines viscous evolution with photoeva poration from the central star 
iClarke. Gendrin & Sotomavoill200il) . In this case, although the 
star spends a relatively long time in the CTTS state (determined 
by the viscous timescale of the outer disc), at late times photoevap- 
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oration becomes important. This is effective only at radii beyond 
about 5 AU and the resulting mass loss from this region deprives 
the inner disc of re-supply. Consequently, the inner disc drains on 
its own (short) viscous time-scale, and the time-scale for the disap- 
pearance of the inner disc is much less than the disc's lifetime to 
date. 

However, as stressed bv [c iarke et al] j2ppi[) (see also 
iMatsuvama. Johnstone & Hartm ann 2003; Ruden 2004), such a 
mechanism can only operate if the photoionizing radiation from 
the central star does not itself decline dramatically as the accretion 
onto the star is shut off. This implies that there must be a significant 
contribution to the photoionizing radiation in T Tauri stars which is 
not powered by accretion if such photoevaporation models are to 
be successful. TTS must not only produce a strong ionizing flux (of 
order lO^'^ ionizing photons per second), but must sustain this flux 
at the late stages of their evolution from CTTS to WTTS. 

Currently, the origin of the photoionizing emission from 
TTS is unclear, and even its magnitude is p oorly constrained 
jGahm et alJll979t llmhoff & ADpe nzeiieilll987t) . The rate of Ly- 
man continuu m pho tons emitted by the solar chromosphere is of 
order 10^**s"^ fflasri et al. 1979; Avre7l99?), and studies of Her- 
big Ae/B e stars indicat e ionizing fluxes of order 10*^-10''^s"i 
I Bouret & Catalal Il998l) . so it seems logical that the value for 
TTS lie between these two. However no firm evidence exists re- 
garding the ionizing con tinuum emitted by TTS. One possible 
emission mecha nism (see Calvet & Gullbrins 1998; Lamzin 199^, 
iGullbring et alii2000) is that it is generated mainly in an accre- 
tion shock where material from the disc impacts the stellar sur- 
face. This "accretion luminosity" has been estimated to be capa- 
ble of providing ionizin g fluxes of order 10^°-10'^^ photon s~ ^ 
faoUenbach et al. 2000 ; Matsuvama et al .''2003; Fo nt et alJl2004b . 
although this does not take into account the fact that the accre- 
tion column itself is highly optical ly thick to ionizing photons 
iAlexander. Clarke & Pringldl2004al> . Further, it is evident that in 
this case the ionizing radiation would decline as the accretion rate 
declined, and phot oevaporation would not be a viable m echanism 
for disc dispersal <Matsuvama et alJl2003l lRudenll2004h . On the 
other hand, other authors have argued that this radiation instead 
derives from a scaled up version of solar-like magnetic activity 
^Lago . Penston & Johnstone 1984; Costa et al. 2000). In this case, 
a decline in accretion rate would not be expected to shut off the 
source of photoionizing radiation, and photoevaporation models for 
disc dispersal would in principle be viable. To date, attempts to dis- 
tinguish between these models on the basis of archival lUE spec- 
tra tCosta et al. 2000; Johns-Krufl, Valenti & Linsky 2000) have 
come to rather different conclusions about the relative contri- 
bution of accretion to the ultraviolet output of T Tauri stars. 
In the X-ray domain, by contrast, it is evident that the main 
energy source is not accretion, since WTTS are, if anything, 
even more luminous in the X-rays than their CTT S counter- 
parts jPamiani etaLlll995tfStezIer & NeuhauseJl200lh . However, 
as long as the central object provides a strong UV radiation field 
X-rays are unlikely to be able to influence disc evolution signif- 
icantly, as the mass-loss rates due to X-ray photoevaporation are 
much sm aller than those driven by the fiduci al UV photoevapora- 
tion rate <AIexander. Clarke & Pringlel2004d) . 

In this paper, we use archival data from HST STIS and from 
the lUE final archive I Valenti, Johns-Krull & Linskv 2000) to re- 
examine the issue of the magnitude and origin of the photoionizing 
emission in T Tauri stars. The structure of the paper is as follows. 
In Section |2| we discuss the problems involved in estimating ion- 
izing fluxes from observations, taking particular note of the large 



uncertainties which arise (primarily due to uncertainties in redden- 
ing). In Section|3|we use an emission measure method to estimate 
the ionizing fluxes from a small sample of CTTS. Following on 
from this (Section |4}, we propose a simple, reddening indepen- 
dent, method of estimating the fraction of the UV power radiated at 
wavelengths short of the Lyman break, namely through the ratio of 
the He II 1640A line to the C IV 155 0A line. Using data from HST 
STIS and from the lUE final archive I Valenti. Johns-Krull & Lins'k^ 
12000.^ we examine how this ratio behaves as TTS evolve. We find 
evidence that the relative contribution of high frequency radiation 
increases, or at least does not decrease, as the accretion rate on 
to the star declines and interpret this as evidence for an important 
chromospheric contribution to the ionizing flux in T Tauri stars. In 
Section|5|we discuss further caveats that apply to our technique, as 
well as some issues regarding source geometry, and summarise our 
conclusions in Section|6| 



2 OBSERVATIONAL PROBLEMS 

There are a number of problems involved in trying to determine 
the nature of the ionizing flux from T Tauri stars. Direct observa- 
tions of their spectra at wavelengths shortward of the Lyman break 
at 912A are not possible as any escaping photons are readily ab- 
sorbed by (abundant) interstellar neutral hydrogen. Observations at 
slightly longer wavelengths (~ 1000-2000A) tend to focus on line 
emission, and consequently have low continuum signal-to-noise. 
This, combined with the large observational uncertainties in de- 
rived stellar parameters (such as distanc e and stellar radius; see, 
for example. iKenvon & Hartmannlll99^ and the uncertainty re- 
garding what is "true" continuum emission and what is unresolved 
line emission, mean that the observed continua are not useful in 
constraining this problem further. It is also impossible to evaluate 
the ionizing flux by studying the B aimer lines of neutral hydrogen, 
which can be due to radiative recombination, as previous studies 
(Valenti. Basri & Johns 1993) have found that these lines are in lo- 
cal thermal equilibrium or, at the very least, that the excitation is 
collisionally dominated. 

Correcting for reddening is also a serious issue in the ultra- 
violet (UV), as extinction is far greater here than in the optical or 
infra-red. Previous observations have resulted in a range of red- 
dening parameters being derived for these objects, with measure- 
ments of Av by different methods commonly differing by a mag- 
nitude or more. (For T Tau observed value s range from Ay = .8 
toavtsova & Lamzirj200l) to Ay = 1.7 l lGullbring et all2000h .) 
When extended to the UV such uncertainties, combined with un- 
certainties regarding the behaviour of the reddening law, result in 
very large uncertainties in the true UV fluxes of T Tauri stars. We 
attempted to obtain independent estimates of the reddening towards 
our HST sources but this proved impossible by conventional meth- 
ods. Interstellar absorption lines such as the S II triplet at 1250A are 
too weak and poorly resolved, and the rapid variation, on time- 
scales of weeks to days, of emission from T Tauri stars means that 
archival optical/UV data cannot be used due to the non-simultaneity 
of the observations. 



3 EMISSION MEASURE ANALYSIS 

In order to make some quantitative estimates of the ionizing fluxes 
emitted by TTS we have first made use of an emission measure 



(EM) calculation. In plasma physics, the differential emission mea- 
sure (DEM) is defined as: 

DEM{T)=n,nH^ (1) 

where rie and uh are the particle number densities of electrons and 
hydrogen respectively, T is temperature and h is length along the 
line-of-sight. Thus the intensity J of a single spectral line due to a 
transition between atomic states i and j can be evaluated from the 
DEM as: 

I = Ah{X) j C{T,\r,,n,)DEM{T)dT (2) 

where \ij is the wavelength of the line, Ab{X) is the abundance of 
the element {X) in question and C is the contribution function, 
which incorporates all the necessary atomic physics. The DEM 
gives an indication of the amount of plasma emitting along the 
line-of-sight in the temperature interval [T, T + dT], and the DEM 
and the necessary atomic data can therefore give intensities for a 
complete spect rum. (For a detai led review of the necessary physics 
see fMason & Monsignori Fossikl994. ) DEMs are derived from ob- 
served spectral data but the process is complex and difficult, requir- 
ing high-quality data. Consequently few DEMs exist for TTS, as 
the spectral data observed are not typically of high enough quality 
to enable a DEM to be derived. 

We have mad e use of the EMs for 5 CTTS derived by 
iBrooks et alJ i200ll) (kindly provided in electronic form by David 
Brooks): BP Tau, RY Tau, RU Lup, GW Ori and CV Cha. The anal- 
ysis of Brooks et al. ( 2001 ) used line fluxes, observed with the lUE 
satellite, to calculate the distribution of the EM. Each spectral line 
results in an "emission measure locus" centred around the forma- 
tion temperature of the line, and so by sampling a series of lines 
over a range of formation temperatures iBrooks et alj i200ll) were 
able to produce an EM distribution as a function of temperature 
for each object studied. The lines used were C II 1335A, the Si IV 
doublet at 1393A, the Civ doublet at 1549 A, Olll] I66OA, Nlll] 
I752A, Sin 1816A, Sim] 1892A and ClIl]1909A. In addition the 
blended Silll/Ol pair at 1303A and the Nv 1242A line, which is 
blended with Lya, were also used. This results in derived EM dis- 
tributions that span the range log T ~ 3.9 -5.6. 

The EM distributions evaluated by IBrooks et alJ i200ll) are 
volume EMs, and so in order to convert them to column EMs we 
have assumed spherical symmetry and scaled the distributions by 
a constant factor 4tv R^/d^ accordingly^. Here R is the stellar ra- 
dius and d the Eart h-star distance, and we have adopted the stel- 
lar parameters from lBrooks & CostaH2003^ . We performed a sim- 
ple polynomial fit to the EM loci to obtain single-valued DEMs as 
functions of temperature. We then used the resulting scaled DEMs 
as inputs to the CHIANTI spectral synthesis code (version 4.01, 
iDere et alJl997l:lYoung et alj2003h . This code makes use of a large 
atomic database to simulate the thermal emission from optically 
thin plasmas. It is useful in our study as it requires few input pa- 
rameters, the most important of which is the EM distribution. In 
this manner, we are able to generate synthetic spectra over a much 
larger wavelength range that can be observed. 

As a test of the reliability of the procedure we have first at- 
tempted to reproduce the observed spectra from the models. In 
order to achieve this we must correct the synthetic line fluxes 
for geometric dilution, again using the stellar parameters from 

^ A volume EM takes the same form as the column EM defined in equation 
m but with the length element dh replaced by a volume element dV. 
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BP Tau: Observed Line Strengths v Synthetic Line Strengths 
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Figure 1. Comparison between observed line strengths and the line 
strengths obtained from the EM m odel for BP Tau. O bserved strengths 
(the extinction-corrected values from lBrooks et alj200lh are shown as error 
bounds, with the synthetic values shown by filled circles. 

lBrooks&Cost3<2003h . We have focused on the strong, unblended 
lines, and found that in general we are able to reproduce the 
observed li ne strengths (the extinction corrected values listed in 
IBrooks eTaL 2001) to within a factor of 2-3. Both the line ratios 
and the absolute strengths are in general well-matched. (The ex- 
ceptions to this are some of the forbidden lines which, as noted 
by Brooks et al. (2001), are extremely density-sensitive and not al- 
ways well matched to the observations.) A comparison of the ob- 
served and calculated line strengths for a typical case (that of BP 
Tau) is shown in FigQ 

We note that the EM distributions are not always well- 
represented by a low-order polynomial, and in some cases the re- 
sulting line-fluxes are very sensitive to the fit. The resultant spec- 
trum also requires that the elemental abundances be specified. 
Again, the line-fluxes are somewhat sensitive to the choice of abun- 
dances, but we are able to reproduce the observed values reason- 
ably well using standard abundance data (the "solar hybrid" abun- 
dances of Fludra & Schmelz 199^. We also note that the predicted 
strength of the C IV 1550A doublet is generally slightly too large, 
which may be an indication of opacity. Lastly we note that the syn- 
thetic spectra completely fail to reproduce the Hell 1640A line, 
with predicted line strengths 2-3 orders of magnitude smaller than 
those o bserved. (This line is not included in FigQas it was not 
used bv lBrookset"an<200ll) in the derivation of the EM and so no 
extinction-corrected value is quoted.) This occurs despite the fact 
that the excitation temperature for this line is very close to that of 
the well-reproduced CiV 1550A line. We attribute this to the fact 
that the observed Hell 1640A line emission is due to radiative re- 
combination, rather than collisional excitation, and will return to 
this point later. 

Whilst the spectra are generally well reproduced this is not 
true for GW Ori, where the predicted absolute line strengths are 
too large by a factor of ~ 100. The error is approximately con- 
stant in all of the line strengths and so it appears to be a systematic 
error, probably due to the scaling of the EM. The quoted radius 
for GW Ori (8.4R0) is somewhat uncertain, but this uncertainty 
is unlikely to be responsible for such a large error. We note that 
Brooks & Costa (200i) also found GW Ori to be anomalous, at- 
tributing their discrepancy to the failure of the assumption of spher- 
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Table 1. Values of the photon fluxes derived for the 5 CTTS studied 
by Brooks et al. (2001) and Brooks & Costa (2003). The fluxes are pre- 
sented in 3 wavebands: ionizing (700-912A), "H2 -dissociating" (912- 
llOOA) and "broadband FUV" (912-2000A). These values were derived 
by adopting the ionization e quilibria from Arnaud & Rothenflua 1 1985) 
and]Arnaud_&_Rayin^^ and "solar hybrid" elemental abundances 

iFluA^^^chmdangggh : different choices can alter these values some- 
what. The extinction parameters Ay adopted in deriving the emission mea- 
sures are included for reference ffrom lBrooks & Costal2003l) . 



Star 


Photon flux in waveband ( 
700-912A($) 912-llOOA 


xl042s-i) 
912-2000A 


Av 


BPTau 


0.68 


1.5 


6.6 


0.5 


RYTau 


2.1 


2.3 


14 


0.55 


RULup 


2.3 


2.6 


16 


0.4 


GWOri 


13 


9.4 


56 


0.8 


CVCha 


78 


71 


440 


1.7 



ical symmetry. In addition we note that GW Ori is a spectroscopic 
binary jMathieu et al. 1991), which may introduce further uncer- 
tainty. As a result we have scaled down the EM for GW Ori by a 
further constant factor of 100 in order to reproduce the observed 
spectrum better. 

As a further check we have attempted to reproduce the values 
of the "total radiated power" derived by Brooks & Costa (2003). 
Again, we can reproduce all these values to within the uncertain- 
ties introduced by the unknown parameters discussed above. Thus 
we are satisfied that the EM distributions and the CHIANTI code 
provide a consistent way of synthesizing the spectra of the chro- 
mospheres of these TTS. 

In this maimer we have used the CHIANTI spectral synthesis 
code to generate synthetic spectra of TTS chromospheres for the 
5 sources mentioned above. As a result, we can use these spec- 
tra to make estimates of the ionizing fluxes they produce. We es- 
timate this by summing the contributions to the synthetic spectra 
over the wavelength range 700-9 12A. Extending the lower limit 
of this range to smaller values makes little difference to the total, 
as the Lyman continuum is dominated by photons at wavelengths 
very close to the Lyman break. For comparison we also evaluate 
the photon fluxes in two wavebands longward of the Lyman break, 
912-llOOA and 912-2000A. The derived values are comparable to 
previous estimates of the FUV flux from TTS (eg. iHerczeg et alj 
[2004). These estimates are subject to the modelling uncertainties 
mentioned above, and more significantly are also subject to the 
much larger uncertainties due to reddening. As a result we consider 
these values to be order-of-magnitude estimates only. However as 
shown in TableQ the resulting ionizing fluxes are high enough to 
meet the demands of disc photoionization models, and in fact may 
be somewhat higher than the previously assumed fiducial value of 
10*^ photon s~^. Thus we find that the Lyman continuum produced 
by TTS may indeed be large enough to drive disc photoionization 
models. However a large ionizing flux alone is not sufficient; we 
also require that the ionizing flux does not decrease as the objects 
evolve. Thus we now seek a diagnostic that will tell us about the 
evolution of the ionizing flux. 



4 A NEW TECHNIQUE 

As noted in Section|3| whilst the synthetic spectra generated by the 
EM analysis are generally accurate, they fail completely to match 
the observed fluxes of the Hell 1640A line. This line is the Hell 



Table 2. Observed He Il/C IV line ratios and derived He Il-ionizing/total 
power ratios for the 3 o bjects where the power ratio is robust. Observed 
line ratios are taken from lValenti et alj <2000t) . 



Star 


He II/C IV 


Power at A < 228A / Total power 


CV Cha 


0.12±0.04 


0.0064 


RULup 


0.15±0.01 


0.0054 


BPTau 


0.50±0.01 


0.0277 



equivalent of Ha, and so may be produced by radiative recombi- 
nation^. This would explain its absence from our synthetic spectra, 
which incorporate only collisional excitation of atoms and ions. 

Bearing this in mind we propose a new technique for deter- 
mining the behaviour of this ionizing flux. If the He II 1640A line 
is due to radiative recombination (there is some debate about this 
- see Section l?^ . then the flux in the He II line should provide an 
indication as to the power of the ionizing radi ation at wavelengths 
< 228A. In addition. Bro oks & Costal ]2003h found that the line 
flux in C IV doublet at 1550A is strongly correlated with the total 
power radiated by TTS atmospheres. As a result, we propose that 
the He II/C IV line ratio should provide an independent indicator of 
the relative flux of photons which are ionizing material around the 
stars in our sample. We h ave investigated co rrelations between the 
observed line ratios (from lValenti et alj200(A and the ratio of He II- 
ionizing power to total power ratios evaluated from the EMs, for 
the 5 objects studied above. Unfortunately the power ratio derived 
from the EM analysis is very sensitive to the slope of the EM dis- 
tributi on in the range log T ~ 4.2^.8. The data from Bro oks et alJ 
teOOl') show some scatter in this region, and as a result some of 
the derived power ratios are very sensitive to the fit. Only BP Tau, 
RU Lup and CV Cha give robust power ratios for various fits. The 
other 2 objects (RY Tau and GW Ori) produce power ratios that 
very uncertain, making the numbers useless in this regard. As seen 
in Table|5|the values for BP Tau, RU Lup and CV Cha do indicate 
that a harder spectrum produces an increased line ratio, supporting 
our hypothesis. 

The advantages of using this line ratio are twofold. Firstly, it 
provides a "normalised" measurement, essentially measuring "ion- 
izing flux/total chromospheric power" and so should be robust 
against variations between sources (which is important, as our sam- 
ple covers a broad range of spectral types and evolutionary states). 
Secondly, the fact that these lines are closely spaced in wavelength, 
separated by only 90A, means that the ratio is relatively insensitive 
to variations in reddening. Indeed we found that adopting extreme 
values of both Ay and Rv (ie. the highest and lowest found in 
the literature) and dereddening observed spectra using the extinc- 
tion curves of lCardelli. Clayton & Mathi3 il989l) caused variations 
of less than 10% in the line ratio: whilst the absolute line fluxes 
are very sensitive to extinction this line ratio is not. In essence we 
propose that this line ratio represents a "spectral hardness ratio", 
related to the fraction of the total chromospheric power radiated at 
short wavelengths. As can be seen from X-ray observations, there 
is no evidence that magnetic activity in TTS is related to disc evo- 
lution, as both CTTS an d WTTS show similar levels of activity 
iFeigelson & Montmerlelll999i) . Thus studying the dependence of 
this line ratio against various evolutionary indicators should give an 
indication as to how the ionizing flux evolves. 



2 The wavelength of the "Lyman break" for He II is 228A. 
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Table 3. Line strengths obtained from HSTSTIS. Values for C IV 1550A are 
tlie sum of tlie two fitted components. 



Hell/CIV Line Raiio v Stellar Age 



Star 


No. of 


Flux (xlO" 


^■^ergs -^cm ^) 


Hell/Civ 




Spectra 


Civ 1550 


Hell 1640 


Ratio 


TTau 


10 


19.94±1.34 


7.01±0.15 


0.35±0.03 


RYTau 


1 


0.73±0.24 


0.17±0.04 


0.23±0.10 


SU Aur 


2 


3.60±0.15 


0.76±0.03 


0.21±0.02 


GW Ori 


2 


2.58±0.52 


0.67±0.02 


0.26±0.05 


COOri 


2 


0.09±0.02 


0.03±0.01 


0.31±0.08 


EZOri 


1 


0.75±0.05 


0.14±0.01 


0.19±0.02 


VI 044 Ori 


1 


1.58±0.09 


0.26±0.02 


0.16±0.02 


P2441 


1 


0.37±0.02 


0.09±0.02 


0.23±0.06 


RYLup 


2 


11.80±0.97 


1.59±0.13 


0.14±0.02 
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Figure 2. Line ratios from STIS plotted against those from lUE, for the 5 
objects observed by both satellites. The solid line represents equality. 



4.1 Observations and Data Reduction 

UV spectra of 9 CTTS were obtained from the Space Telescope 
Imaging Spectrograph (STIS) on board the Hubble Space Tele- 
scope (HST), taken from the H ST public archive (see Table|3|for the 
names of these objects). Initial processing (flatfielding, cosmic ray 
removal etc.) was perform ed by the HST STIS reduction pipeline 
iTravisano & R ichon* 1997) and so we were provided with flux- 
calibrated 2-dimensional spectra. These spectra were added across 
the spatial extent of the targets and a sky subtraction was performed 
using an offset of 2'. In cases where more than one spectrum was 
observed (see Table |3}, the spectra were co-added to improve the 
S/N ratio. The fluxes for the He II and C IV lines were then extracted 
using the line-fitting routines from the FIGARO package. The He II 
I640A was fitted using a single Gaussian profile, whereas the C IV 
line was fitted by 2 separate Gaussian components corresponding 
to the 1548A and 155lA peaks. The results of these line fits are 
also shown in Table|3| 

The sample of 9 objects observed by STIS is small, and also 
covers a small spread in evolutionary properties, as all of the 9 ob- 
jects are CTTS and have similar ages. Further, the objects cover a 
broad range of spectral types (F8-K4) and so this sample is far from 
ideal. In order to expand the sample we turned to dat a from the In- 
ternat ional Ultraviolet Explorer (lUE) final archive iValenti et alj 
l200(t) . which observed a total of 50 TTS, both classical and weak- 
lined. In some cases the signal-to-noise ratio is very poor, so only 



Figure 3. Plot of line ratio versus stellar age. (Ages taken from 
IPalla & StahleJ2002l l 



Hell/CIV Line Ratio v Accretion Rate 



dU/dt / Mg^r 

Figure 4. Plot of line ratio versus accretion rate: objects "evolve" to the 
right. (M values from .Gullbring et al. 19 98) 

33 of these objects have reliable line fluxes for both C IV 1550A and 
Hell 1640A. The lower spectral resolution of lUE (compared to 
that of STIS) intro duces some p r oblem s with line-blending espe- 
cially, as noted by ' Valenti et alJ llOO& i. in the Hell 1640A fine. 
However, plotting the Hell/Civ line ratio from STIS versus that 
from lUE, as shown in Fig|2| seems to indicate that the line blends 
provide a systematic offset to this ratio. If the variations were due to 
variability we would not expect to see such a systematic effect, and 
so it seems that whilst the absolute values of the line ratio differ be- 
tween the two samples the appearance of any evolutionary trends 
should be not be affected by this systematic effect. Consequently 
we use the ratios obtained from the lUE archive throughout the re- 
sults section, as they are as valid as those from STIS and incorporate 
a much larger sample. 



4.2 Results 



Figures 13 15 1 show the He II/C IV line ratio (from the data in 
f/alenti et al. 2OO0) plotted against several evolutionary indicators 
(where available). While good values of the line ratio are available 
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Hell/CIV Line Ratio v Ha Equivalent '.Vidtti 



• CTTS 

o vms 



EW(Hei) / S 



Figure 5. Plot of line ratio versus Ha e quivalent width: obiec ts "evolve" to 
the right. (Equivalent widths taken fromlCoh en & Kuhll979t ^ SU Aur is a 
somewhat unusual object i Unruh et al. 2004), but it does posess a disc and 
so we have classified it as a CTTS despite its weak Ha equivalent width of 
3A. 



for 33 objects, in some cases the evolutionary indicators are less 
well-known. Thus not all of the objects appear on all of the plots. 

Fig. |3| shows the line ratio plotted against stellar age, with 
values for the age s, derived from pre-ma in-sequence evolutionary 
tracks, taken from lPalla & Stahie3i2002h . No clear trend is visible, 
although there is a slight trend for a higher ratio in older systems. It 
may also be significant that the highest values of the line ratio are 
both for WTTS. 

FigElshows the line ratio plotted against mass accretion rate. 
Mass accretion rate is difficult to measure, and only a small sam- 
ple of objects have known acc retion rate s. We have used the ac- 
cretion rates derived by Gu Ubring et alj fl998). derived by fit- 
ting accretion-shock models to observed spe ctra. Other similarly- 
sized samples of mass accretion rates exist faartiganetai]|l995t 
Ijohns-Kru U et al.'2000Y but the systematic errors between the dif- 
ferent methods used to derive the accretion rates mean that the data 
cannot be combined. Again no strong trend is visible, but if any- 
thing the line ratio appears to increase as the objects evolve. 

FigOshow s the line ratio p lotted against Ha equivalent width, 
as measured by Cohen & Kuhi ( 1979). The Ha line arises due to 
accretion onto the central object, and thus decreasing Ha e mission 
is thought to be a good "disc clock" jKitamura et all2002h . Again, 
no clear trend is visible, but we note again that the highest values 
tend to be for the WTTS. 

In all cases the line ratio appears to increase, or at least not to 
decrease, as the objects evolve. This appears to imply that the ion- 
izing flux is not powered by mass accretion, as the fraction of the 
total flux emitted at short wavelengths, and consequently the ob- 
served He II/C IV line ratio , would be expected to decrease dramat- 
ically in this case ( Matsuv ama et all2 003'). Accretion shock models 
iCalvet & Gullbrina„1998f) predict that the emission will decrease 
and peak at longer wavelengths as individual objects evolve, al- 
though an increased mass-to-radius ratio results in an intrinsically 
harder emitted spectrum. However no correlation of any kind is ob- 
served between M/R and the Hell/Civ line ratio and so we do 
not consider this to be a significant effect; both high and low line 
ratios are obtained from objects with both high and low M/ R. It is 
also significant that the highest values of the line ratio are obtained 



for the WTTS in the lUE sample, as these are expected to be disc- 
less objects and so disc accretion cannot be happening here. Thus it 
seems likely that the production of ionizing photons around TTS is 
a stellar phenomenon, which we attribute to the chromosphere, and 
is independent of disc accretion. As mentioned in Section|4| there 
is no evidence that magnetic activity in TTS is age-dependent, and 
so a large ionizing flux fraction should imply a strong absolute ion- 
izing flux. Thus it seems likely that TTS can produce an ionizing 
flux sufficient to power disc photoevaporation models: the ionizing 
flux is strong, and shows no significant decrease with time. 



5 DISCUSSION 

5.1 Emission Measures 

The use of the emission measure analysis is obviously crude, and 
contains many uncertainties. The dominant uncertainty in the abso- 
lute values of $ derived is the reddening uncertainty. The manner in 
which the EMs are scaled also results in a form that is weakly sen- 
sitive to the choice of stellar parameters. However this uncertainty 
is small compared to that introduced by reddening. 

As regards the "hardness" of the spectrum, specifically the rel- 
ative power radiated longward and shortward of the Lyman break, 
this is also not especially well constrained. As mentioned in Section 
|4|above, the ionizing/total power ratios derived are rather sensitive 
to the EM fits, and so the correlation noted in table|2|is somewhat 
dubious. We note however that the absolute values of "J> are not es- 
pecially sensitive to the fit, with the fitting uncertainty remaining 
much less significant that that due to reddening. 

However, we are able to reproduce both the observed line 
strengths, both in absolute and relative terms, to a good degree of 
accuracy. Consequently we are satisfied that the EMs are consistent 
with the observed data, although we note that they are not always 
strongly constrained. The absolute values of $ derived in Section|3| 
are considered to be accurate only to around an order of magnitude, 
but they still provide a valuable new constraint. 



5.2 He II Line 

There are several obvious problems with the use of this line to 
infer the behaviour of the ionizing flux from TTS. Firstly, and 
most importantly, the Hell ionization energy is 54.4eV, and so 
recombination lines can only be used as an indicator of the flux 
of photons with energies greater than this (ie. with wavelengths 
< 228A). Consequently this measure omits ionizing photons with 
wavelengths between 228 A and 912A, and photons in this wave- 
length range dominate the ionization of neutral hydrogen. However 
it is difficult to envisage an emission process which could produce 
photons selectively in this wavelength region without producing 
higher-energy photons also, and so whilst the He II recombination 
line may not be ideal it should provide a reasonable indication of 
the flux ionizing neutral hydrogen. 

Secondly there is the issue of the He II line itself, and whether 
it arises from radiative recombination (RR) or collisional excitation 
(CE). This line has been observ ed in the so lar coron a, and has bee n 
explained both in terms of RR ( Zirin"l97^ and CE ('Jordai ll973) . 
[Linskv et al. ( 1995, 1998) have studied the He II 1640A line in the 
chromosphere of Capella and found evidence for both mechanisms, 
but used high-resolution spectroscopy to show that the line profiles 
are ch aracteristic o f RR being the dominant factor. More qualita- 
tivelv. lzirii] il975h found that his RR model works extremely well 
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in active regions of tiie solar chromosphere, but less well in the 
quiet corona; given that T Tauri stars are characterised by their high 
levels of chromospheric activity it seems that RR will play a sig- 
nificant role. Our synthetic spectra dramatically under-estimate the 
strength of the He II 1640A line whilst reproducing the other lines 
and the continuum well, and whilst this lends further support to the 
RR theory we are unable to say with certainty whether the emis- 
sion is CE or RR in origin. However, it is worth bearing in mind 
that the densities and temperatures of plasma required to produce 
significant CE emission from He II are such that we would expect 
such a plasma to radiate significantly at wavelengths short of the 
Lyman break through free-free and bound-free transitions. As a re- 
sult, even if the emission is not RR, the presence of significant He II 
1640 A emission should be a reliable indicator of a significant ion- 
izing flux. The 1640A line is strong in almost all of the sources, 
both young and old, classical and weak-lined, and so it does seem 
that a significant ionizing flux is produced by all of these objects. 

One should also note the fact that the ionization energy for 
He II is 4 times greater than that for H I. As a result the strength of 
the He II 1640 A line may be due to ionization by X-rays rather than 
softer ultraviolet photons. Indeed, there is some evidence for a cor- 
relation between the He II 1640A flux and X-ray power ( Costa et al. 
[2000). However this may merely indicate that both are good tracers 
of magnetic activity, rather than suggesting that the X-rays stimu- 
late the He II 1640A emission. We are not able to make any distinc- 
tion in this regard, but the conclusion that He II 1640A line emis- 
sion traces the production of energetic photons seems secure. 



5.3 C IV Line 

There is also some uncertainty in the use of the Civ 1550A line 
as a normalising factor, as the relationship betwe en C IV flux and 
total r adiated power is purely an empirical one iSrooks & Costal 
i2003h ). As mentioned above, the total radiated power is evalu- 
ated from an emission measure integral covering the UV and X-ray 
regimes. They find that the X-ray contribution to this flux is "less 
than 10%" of the total, but they do not state whether the ionizing 
or no n-ionizing UV flux dominates the total. However a previous 
paper iCosta et al . 2000) indicates that the contributions are similar 
and that, if anything, the non-ionizing UV flux is more significant. 
Our emission measure analysis supports this, as the ionizing flux is 
not the dominant contribution to the total for any of the sources we 

consider (see Table0. 

It should also be noted that iBrooks & Costal j2003) assume 
the C IV emission to be optically thin. If this is the case the rela- 
tive intensities of the two components of the C IV line should differ 
by a factor of 2 (Ardilaetal. 2002). We observe some evidence 
for opacity in 3 of the 9 objects observed with STIS (RY Tau, SU 
Aur and EZ Ori show ratios in the range 1.5-1.8). We also noted 
possible evidence for opacity when comparing our synthetic spec- 
tra to the observed data (see Section|3j, although this may not be 
significant. Therefore there is an unquantifiable degree of uncer- 
tainty here, as the resolution of the lUE data does not permit such 
analysis to be carried out. We also note that a fraction of the C IV 
line flux may be due to accretion, and so the slightly higher line ra- 
tios observed in the WTTS may not indicate a higher ionizing flux. 
However the C IV 1550A line is one of the main coolants in T Tauri 
atmospheres/stellar chromospheres and so it is at least qualitatively 
reasonable to expect the relationship b etween C IV flux and total 
radiated power to hold. Moreover, Broo ks & Costal |2003) found a 
strong correlation over nearly 2 orders of magnitude in power, and 



so we feel that the uncertainties due to the C IV line are less of a 
concern than those due to the He II line discussed in Section l5!2l 



5.4 Geometric issues 

In addition to merely identifying the presence of a strong ionizing 
flux, there is also the problem of getting the photons to the disc. As 
we have previously noted, the accretion col umns are extremely op - 
tically thick to Lyman continuum photons ( Alexander et alj2004al) . 
and so it seems likely that any ionizing photons emitted by the chro- 
mosphere at low latitude will likely be absorbed by the columns. In 
addition, any (bi-)polar outflow may also be optically thick to ion- 
izing photons (Shane et al. 2002). However in order to influence 
disc evolution it is only necessary to provide a strong ionizing flux 
to the disc at late times, and it is not clear if any outflow will per- 
sist to this stage of the evolution. We also note that our diagnostic 
relies on a He II recombination line rather than a H I recombina- 
tion line, so we may see He II emission from regions where Lyman 
continuum emission is suppressed by absorption. 

However there are also positives for the photoevaporation 
model. Chromospheric emission at high latitude and/or from sig- 
nificantly above the photospheric surface will likely be unaffected 
by the accretion columns and, as mentioned above, outflows are 
much less signific ant at late times. Also, wh ilst non-axisymmetric 
accretion columns ^ Mahdavi & Kenvonll998,) will still absorb ion- 
izing photons from the hotspot at their base they will permit much 
of the (roughly symmetric) chromospheric emission to escape un- 
harmed. Thus whilst the TTS geometry is not entirely favourable, 
it seems likely that most reasonable geometries will allow a sig- 
nificant fraction of the chromospheric ionizing flux to escape and 
influence the disc at late times. 



6 SUMMARY 

By analysing data from UV studies of TTS we have placed new 
constraints on the ionizing flux emanating from them. An emission 
measure analysis using existing data has found rates of ionizing 
photons from the chromosphere in the range ~ 10*^-10^* pho- 
tons s^^forSCTTS with ages of around lO^yr. It should be noted, 
however, that these values are subject to large observational uncer- 
tainties, in particular due to reddening. We also note that our de- 
rived ionizing fluxes are intermediate between the solar value and 
the values typical of Herbig Ae/Be stars, and so seem qualitatively 
reasonable. 

We have proposed that the He II/C IV line ratio can be used as 
a reddening-independent diagnostic of the ionizing flux from TTS, 
essentially measuring the hardness of the spectrum. We have anal- 
ysed the behaviour of this line ratio in a larger sample of TTS ob- 
served by lUE, and find no evidence of any correlation with stellar 
or disc evolution. Thus we conclude that the observed emission 
is likely due to the central object, presumably the stellar chromo- 
sphere, and is not tied to disc accretion. This supports the hypoth- 
esis that the chromospheres of TTS can provide a strong ionizing 
flux. Our method is somewhat crude and the results are by no means 
conclusive, but we have begun to place some constraints on a prob- 
lem which was previously unbounded by observations. 

These chromospheric Lyman continu a appear to be sufficient 
to dri ve models of disc photoevapo ration iHollenbach et alj|l994t 
Clark e et alJl200ll:lFont et alJl2004 . However additional observa- 
tions in the UV are clearly still needed to resolve this issue, as the 
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current sample is small and the data are subject to large uncertain- 
ties. In particular, studies of WTTS will be important in further 
constraining this problem, as all of our derived values of the ioniz- 
ing flux are for CTTS. The presence of a strong Lyman continuum 
around WTTS would provide further support for the hypothesis that 
the evolution of TTS discs is strongly affected by photoevaporation. 



ACKNOWLEDGEMENTS 

We thank Bob Carswell for useful discussions. We are grateful to 
David Brooks for providing his emission measure data in elec- 
tronic form. RDA acknowledges the support of a PPARC PhD 
studentship. CJC gratefully acknowledges support from the Lev- 
erhulme Trust in the form of a Philip Leverhulme Prize. RDA and 
CJC thank UC Santa Cruz for hospitality during the early stages of 
this work. JEP thanks STScI for continued hospitality vmder their 
Visitors' Program. CHIANTI is a collaborative project involving the 
NRL (USA), RAL (UK), and the Universities of Florence (Italy) 
and Cambridge (UK). Based in part on observations made with 
the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Science Institute. STScI is operated 
by the Association of Universities for Research in Astronomy, Inc. 
under NASA contract NAS 5-26555. We thank an anonymous ref- 
eree for several comments which improved the clarity of the paper. 



REFERENCES 

Adams, F.C., Hollenbach, D., Laughlin, G., Gorti, U., 2004, ApJ, 

611, 360 

Alexander, R.D., Clarke C.J., Pringle, I.E., 2004a, MNRAS, 348, 
879 

Alexander, R.D., Clarke C.J., Pringle, I.E., 2004b, MNRAS, 354, 
71 

Arnaud, M., Raymond, J.C., 1992, ApJ, 398, 39 

Arnaud, M., Rothenflug, R., 1985, A&ASS, 60, 425 

Ardila, D.R., Basri, G., Walter, P.M., Valenti, J.A., Johns-KruU, 

CM., 2002, ApJ, 566, 1100 
Armitage, RJ., Clarke, C.J., Tout, C.A., 1999, MNRAS, 304, 425 
Ayres, TP, 1997, J. Geophys. Res., 102, 1641 
Basri, G.S., Linsky, J.L., Bartoe, J.-D.F., Brueckner, G., Van 

Hoosier, M.E., 1979, ApJ, 230, 924 
Beckwith, S.V.W., Sargent, A.I., Chini, R.S., Gusten, R., 1990, 

AJ, 99, 924 

Bouret, J.-C, Catala, C, 1998, A&A, 340, 163 

Brooks, D.H., Costa, V.M., 2003, MNRAS, 339, 467 

Brooks, D.H., Costa, V.M., Lago, M.T.V.T., Lanzafame, A.C., 

2001, MNRAS, 327, 177 
Calvet, N., Gullbring, E., 1998, ApJ, 509, 802 
Cardelh, J.A., Clayton, G.C., Mathis, J.S., 1989, ApJ, 345, 245 
Clarke, C.J., Gendrin, A., Sotomayor, M., 2001, MNRAS, 328, 

485 

Cohen, M., Kuhi, L.V, 1979, ApJS, 41, 743 
Costa, V.M., Lago, M.T.V.T., Norci, L., Meurs, E.J.A., 2000, 
A&A, 354, 621 

Damiani, R, Micela, G., Sciortino, S., Hamden, F.R., Jr., 1995, 
ApJ, 446, 331 

Dere, K.P, Landi, E., Mason, H.E., Monsignori Fossi, B.C., 

Young, RR., 1997, A&ASS, 125, 149 
Feigelson, E.D., Montmerle, T., 1999, ARA&A, 37, 363 
Fludra, A., Schmelz, J. T., 1999, A&A, 348, 286 



Font, A.S., McCarthy, I.G., Johnstone, D., Ballantyne, D.R., 2004, 
ApJ, 607, 890 

Gahm, G.F., Fredga, K., Liseau, R., Dravins, D., 1979, A&A, 73, 
L4 

Gullbring, E., Hartmann, L., Briceno, C, Calvet, N., 1998, ApJ, 

492, 323 

Gullbring, E., Calvet, N., Muzerolle, J., Hartmann, L., 2000, ApJ, 
544, 927 

Haisch, K.E., Lada, E.A., Lada, C.J., 2001, ApJ, 553, L153 
Hartigan, P., Hartmann, L., Kenyon, S.J., Strom, S.E., Skrutskie, 

M.R, 1990, ApJ, 354, L25 
Hartigan, R, Edwards, S., Ghandour, L., 1995, ApJ, 452, 736 
Hartmann, L., Calvet, N., Gullbring, E., D'AIessio, P., 1998, ApJ, 

495, 385 

Herczeg, G.J., Wood, B.E., Linsky, J.L., Valenti, J.A., Johns- 
KruU, CM., 2004, ApJ, 607, 369 

Hollenbach, D., Johnstone, D., Lizano, S., Shu, F., 1994, ApJ, 
428, 654 

Hollenbach, D.J., Yorke, H.W., Johnstone, D., 2000, in Protostars 
& Planets IV, eds. V. Mannings, A.P. Boss, S.S. Russell, Univer- 
sity of Arizona Press, p401 
Imhoff, C.L., Appenzeller, I., 1987, in Astrophysics and Space 
Science Library, Vol. 129, Exploring the Universe with the WE 
Satellite, ed. Y. Kondo (Dordrecht: Reidel), 295 
Johns-KruU, CM., Valenti, J. A., Linsky, J.L., 2000, ApJ, 539, 815 
Jordan, C, 1975, MNRAS, 170, 429 
Kenyon, S.J., Hartmann, L., 1995, ApJs, 101, 117 
Kitamura, Y, Momose, M., Yokogawa, S., Kawabe, R., Tamura, 
M., Ida, S., 2002, ApJ, 581, 357 
Kravtsova, A.S., Lamzin, S.A., 2002, Astronomy Letters, 28, 835 
Lago M.T.V.T, Penston M.V., Johnstone R., 1984, in Roffe E., 
Battrich B., eds, Proc. Fourth European lUE Conference, ESA 
SP-218, p.233 
Lamzin, S.A., 1998, Astronomy Reports, 42, 322 
Linsky, J.L., Wood, B.E., Judge, P., Brown, A., Andrulis, C, 

Ayres, T.R., 1995, ApJ, 442, 381 
Linsky, J.L., Wood, B.E., Brown, A., Osten, R.A., 1998, ApJ, 492, 
767 

Mahdavi, A., Kenyon, S.J., 1998, ApJ, 497, 342 

Mason, H.E., Monsignori Fossi, B.C., 1994, A&AR, 6, 123 

Mannings, V, Sargent, A.I., 1997, ApJ, 490, 792 

MaUiieu, R.D., Adams, EC, LaUiam, D.W., 1991, AJ, 101, 2184 

Matsuyama, I., Johnstone, D., Hartmann, L., 2003, ApJ, 582, 893 

Palla, R, Stabler, S.W., 2002, ApJ, 581, 1194 

Richling S., Yorke, H.W., 1997, A&A, 327, 317 

Ruden, S.P, 2004, ApJ, 605, 880 

Shang, H., Glassgold, A.E., Shu, RH., Lizano, S., 2002, ApJ, 564, 
853 

Skrutskie, M. R, Dutkevitch, D., Strom, S.E., Edwards, S., Strom, 

KM., Shure, M.A., 1990, AJ, 99, 1187 
Stezler, B., Neuhauser, R., 2001, A&A, 377, 538 
Strom, K.M., Strom, S.E., Edwards, S., Cabrit, S., Skrutskie, M.R, 

1989, AJ, 97, 1451 
Travisano, J.J., Richon, J.G., 1997, ADASS VI, A.S.R Conference 

Series, Vol. 125, G. Hunt and H.E. Payne, eds., 286 
Unruh, YC. et al., 2004, MNRAS, 348, 1301 
Valenti, J.A., Basri, G., Johns, CM., 1993, AJ, 106, 2024 
Valenti, J.A., Johns-KruU, CM., Linsky, J.L., 2000, ApJS, 129, 

399 

WyaU M.C, Dent W.R.R, Greaves J.S., 2003, MNRAS, 342, 876 
Young, P.R., Del Zanna, G., Landi, E., Dere, K.P., Mason, H.E., 
Landini, M., 2003, ApJS, 144, 135 



Photoionization rates around T Tauri Stars 

Zirin, H., 1975, ApJ, 199, 63 



